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摘 E: 目前 ,冰川 分 布 与 变化 的 影响 机 理 还 存在 不 确定 机 制 。 为 探究 地 形 因 子 对 冰川 分 布 格局 的 影响 ,本 文 以 
第 6 版 Randolph 冰川 编目 数据 和 NASADEM 为 数据 源 ,运用 均值 变 点 法 .分 布 指 数 和 地 理 探测 需 等 方法 ,对 青藏 高 


原 冰川 在 地 形 起 伏 度 ,海拔 ,坡度 和 坡 向 等 地 形 梯度 上 的 分 布 格局 及 其 与 地 形 因子 之 间 的 响应 关系 进行 分 析 。 结 


拔 上 呈现 优势 分 布 。(2) 地 形 对 冰川 发 育 影响 较 大 ,各 地 


果 表 明 :(1) 青藏 高 原 冰川 对 坡 向 \ 坡 度 的 适宜 性 较 广 ,对 地 形 起 伏 度 、 海 拔 的 选择 性 较 强 , 且 在 高 地 形 起 伏 度 、 高 海 


因子 对 冰川 空间 分 布 的 影响 显著 不 同 。 海 拔 和 地 形 起 伏 


度 是 控制 冰川 分 布 的 主要 因素 , 且 交 互 探测 结果 表明 二 者 共同 作用 对 冰川 空间 分 异性 影响 最 大 。(3 ) 在 海拔 和 地 形 
起 伏 度 梯 度 上 ,喜马拉雅 山脉 冰川 优势 分 布 范围 最 广 , 其 次 是 喀 喇 昆仑 山脉 ,其余 10 座 山脉 冰川 的 优势 分 布 均 呈 近 


似 正 态 分 布 ,但 不 同 山脉 形态 各 异 。 


关键 词 : 冰川 分 布 ; 地 形 梯度 ; 分 布 指数 ; 地 理 探测 需 ; 


冰川 是 冰冻 圈 的 重要 组 成 部 分 ,与 地 表 能 量 平 
f .海平 面 上 升 `. 区 域 气候 和 水 循环 等 有 着 密切 联 
AU ,冰川 的 存在 有 助 于 稳定 气候 ,改善 水 资源 , 保 
护 生态 系统 和 地 质 形 态 等 。 青 藏 高 原 是 世界 海拔 
最 高 面积 最 大 的 内 陆 高 原 , 其 独特 的 地 势 高 .气温 
低 且 异 夜 温差 大 的 条 件 非常 有 利于 冰川 发 育 ,作为 
中 低 纬度 现代 冰川 的 作用 中 心 ,也 是 亚洲 乃至 北 半 
球 环境 变化 的 调控 器 。 

冰川 的 形成 与 发 育 取决 于 地 形 条 件 和 降水 A 
温 等 气候 条 件 的 组 合 “。 冰 川 的 分 布 与 变化 目前 
还 存在 不 确定 的 机 制 。 大 多 数 人 研究 将 亚洲 高 山地 
区 的 冰川 物质 平衡 归 因 于 气候 变化 一 ,而 忽视 了 其 
复杂 地 形 的 影响 ”。 然 而 ,由 于 区 域 气候 与 地 形 的 


青藏 高 原 


作用 ; 张 鲜 稚 等 ”利用 第 二 次 冰川 编目 数据 ,研究 
了 中 国 西部 各 个 山脉 冰川 面积 在 海拔 梯度 上 的 分 
布 特征 、 冰 川 平衡 线 高 度 场 的 分 布 规律 和 积累 区 比 
率 的 分 布 特征 ; 周 远 刚 等 研究 发 现 ,高 程 坡 度 和 
坡 向 对 冰川 和 积 雪 面积 变化 的 影响 较 显著 ; 黄 丹 妮 
等 ”对 东 帕 米尔 高 原 冰 川 运动 特征 进行 分 析 ,发 现 
冰川 运动 速度 与 冰川 规模 .坡度 BE Ti], A o PH 
关 ; 张 廉 卿 等 “研究 念 青 唐古拉 山西 段 小 冰期 以 来 
冰川 变化 ,发 现 冰川 的 面积 、 坡 度 和 高 程 对 冰川 面 
让 变 化 的 解释 率 为 71%; 李 林 风 等 “运用 地 理 探测 
带 对 石 羊 河流 域 冰 川 变化 与 地 形 因 子 之 间 的 关系 
进行 研究 ,得 到 海拔 是 冰川 变化 空间 分 蜡 最 重要 的 
控制 因素 。 综 上 所 述 ,前 人 的 研究 考虑 的 地 形 因子 


差异 ,冰川 分 布 与 变化 具有 显著 的 区 域 差异 性 。 深 
入 分 析 这 种 复杂 地 形 的 影响 对 于 理解 冰川 分 布 与 
变化 的 空间 分 异 至 关 重 要 ”" 。 张 正 勇 等 "通过 构 
建 冰川 -地 形 -气候 遥感 反 演 模型 ,对 新 疆 玛 纳 斯 河 
流域 的 冰川 分 布 格局 对 地 理 因子 响应 机 制 进行 研 
究 ,发 现 地 形 因子 对 冰川 的 发 育 和 积累 具有 决定 性 
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多 为 海拔 坡度 和 坡 向 。 山 体 海拔 越 高 ,气温 越 低 ， 
冷 储 越 大 ,加 之 其 对 水 汽 的 拦截 ,更 有 利于 冰川 的 
发 育 ,往往 能 发 育 规模 较 大 的 冰川 ”…" ;缓坡 有 利 
于 冰川 的 累积 和 稳定 ,从 而 可 形成 面积 更 大 的 单条 
冰川 "”。 在 其 他 条 件 相 同 的 情况 下 , 阴 坡 接受 大 
阳 辐 射 热量 少 ,冰川 消融 速率 较 低 ,更 适 于 冰川 发 


基金 项 目 : 国家 自然 科学 基金 项 目 (41671072) ;商丘 师范 学 院 博 士 人 才 引 进项 目 50018401); 商丘 市 科技 创新 领军 人 才 项 目 SQRC202212005); 


商丘 师范 学 院 黄河 故道 生态 保护 和 高 质量 发 展 项 目 (2021KYFZ06) 


作者 简介 : 许 宁 (1987-), 男 ,博士 ,讲师 ,主要 从 事 生态 环境 遥感 与 GIS 应 用 研究 . E-mail: sqnugis@126.com 


通讯 作者 : 李 治 国 . E-mail: lizhiguo999999@163.com 


230 一 239 页 


http: //azr.xjegi.com 


HU, miu JÉ E X RE (Relief Degree of Land Sur- 
face, RDLS) 作 为 地 貌 类 型 划分 的 重要 依据 ,对 评价 
区 域 人 居 环 境 适 宜 性 与 资源 环境 承载 力也 极其 关 
键 “" ,目前 , 鲜 有 地 形 起 伏 度 与 冰川 分 布 与 变化 关 
系 的 研究 。 此 外 ,之 前 的 研究 大 多 利用 常规 统计 回 
归 方 法 研究 冰川 分 布 变化 与 地 形 因子 之 间 关 系 , 但 
该 方法 存在 因数 据 原因 导致 解释 力 不 足 的 劣势 。 
而 擅 于 探测 空间 分 异性 、 揭 示 驱 动力 的 新 统计 方 
法 一 一 地 理 探 测 融 ,同时 具有 可 以 探测 不 同 数 据 类 
型 以 及 探测 两 因子 交互 作用 于 因 变 量 的 优势 ”。 

本 文 以 青藏 高 原 为 研究 区 ,以 第 6 版 Randolph 
冰川 编目 资料 和 NASADEM 为 数据 源 , 除 考 虑 海拔 、 
坡度 和 坡 向 等 地 形 因 子 外 ,引入 地 形 起 伏 度 因子 ， 
并 借助 分 布 指数 和 地 理 探 测 器 技术 对 青藏 高 原 冰 
川 在 地 形 起 伏 度 海拔、 坡度 . 坡 向 等 地 形 因子 梯度 
上 的 分 布 格局 及 其 对 地 形 因子 的 啊 应 机 制 进行 探 
究 ,以 期 为 区 域 水 文 研 究 提供 数据 支撑 和 理论 
参考 。 


1 研究 区 与 数据 源 


1.1 研究 区 概况 
青藏 高 原 边 界 依据 Zhang 等 "青藏 高 原 界线 ， 
南 自 喜马拉雅 山脉 南 药 , 北 至 帕 米 尔 高 原 北 缘 \ 西 
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米尔 高 原 西 缘 , 东 抵 祁 连 山 .横断 山 等 山脉 东 缘 ,总 
面积 为 308.34x10'km*( 图 1)。 青 藏 高 原平 均 海拔 约 
4320 m, 地 势 西高 东 低 ,边缘 山脉 起 伏 程度 较 大 ,内 
部 柴 达 木 盆地 与 芜 塘 高 原 起 伏 程 度 较 小 。 青 藏 高 
原 是 中 低 纬 度 现代 冰川 的 作用 中 心 ,被 誉 为 “亚洲 
水 塔 ”, 且 青藏 高 原 的 冰川 大 都 处 于 高 海拔 低 纬 度 
地 区 ,其 对 于 气候 变化 尤为 敏感 "。 
1.2 数据 来 源 

本 文 使 用 的 冰川 数据 为 第 6 版 Randolph 冰川 编 
H (RGI 6.0) 数 据 ” ,来源 于 美国 国家 雪 冰 数据 中 心 
(http:/www.glims.ore/RGI)。 该 数据 与 之 前 版 本 相 
比 ,提升 了 有 了 明确 时 间 信 息 的 冰川 占 比 ,吸纳 了 区 
域 性 冰川 编目 , 细 化 了 冰川 区 范围 的 界定 。 根 据 
RGI 6.0 的 数据 ,青藏 高 原 地 区 冰川 共 77325 条 ,总 
面积 8.3x10' km。 本 研究 主要 依据 第 6 版 Randolph 
冰川 编目 数据 和 第 二 次 中 国 冰川 编目 ,将 青藏 高 
原 地 区 所 有 冰川 划分 为 12 个 山地 冰川 ,其 中 , 喀 喇 
昆仑 冰川 的 划分 还 参考 了 Mole 等 ”的 研究 。NA- 
SADEM 数据 来 源 于 美国 国家 航空 航天 局 (https:// 
earthdata.nasa.gov/) ,该 数据 是 SRTM 数据 的 升级 版 ， 
空间 分 辩 率 为 30 m7", A&F NASADEM ,利用 Arc- 
GIS 软件 提取 了 青藏 高 原 地 形 起 伏 度 ,海拔 坡度 、 
坡 向 4 个 地 形 因子 。 
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图 1 青藏 高 原 相 对 高 差分 布 


Fig. 1 Distribution map of relative elevation difference of Tibetan Plateau 
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2 研究 方法 


2.1 地 形 起 伏 度 的 计算 

地 形 起 伏 度 "表征 高 程 起 伏 状 况 ,其 计算 存 
在 尺度 依赖 性 ,关键 是 确定 最 佳 分 析 窗 口 。 地 形 起 
伏 度 的 计算 公式 如 下 : 

RDLS = ALT/1000 +{ [Max(H) - Min(H)] x 
[1 — P(A)/A] /500 
式 中 :RDLS 为 地 形 起 伏 度 ;4L7 Max GT) „Min (H) 、 
P(4) 分 别 为 以 某 一 栅 格 单元 为 中 心 的 最 佳 分 析 窗 
口 下 平均 海拔 (四 、 最 高 海拔 (m)、 最 低 海拔 (m) 和 平地 
(相对 高 差 <30 m HIKER) R (km?) A 为 区 域 总 面积 
(kn?) , 

本 文采 用 均值 变 点 法 来 确定 最 佳 分 析 窗 口 “”。 
首先 ,利用 ArcGIS 软件 统计 分 析 窗 口 (2x2,3x3,…， 
121x121) 下 的 平均 高 差 , 再 依次 计算 单位 地 势 度 7， 
公式 如 下 : 


(1) 


T,=t,/s,. (522,3, +++, 120) (2) 
式 中 :7 为 第 个 分 析 窗 口 下 单位 地 势 度 ;t 为 平均 
高 差 ;si 为 分 析 窗 口 面 积 。 对 数列 7T 取 对 数 (1n7) ,得 
到 非 线 性 数列 样本 读 (i=2,3,…,120); 并 将 乱 分 为 
HWA PUBS, BIX, X5, XA Xie Xoo °° ,Xiw, 分 别 计算 
总 样本 的 离 差 平方 (5) 与 前 后 两 段 样本 的 离 差 平 方 和 
(5;)。 二 者 之 差 的 最 大 值 (5-5;) 对 应 的 分 析 窗 口 即 为 
最 佳 分 析 窗 口 。 基 于 NASADEM 数据 运用 均值 变 点 
法 计算 青藏 高 原 地 形 起 伏 度 最 佳 分 析 窗 口 为 36x36 
个 像 元 的 矩形 邻 域 ,对 应 面积 约 为 1.17 km?( 12) 
2.2 分 布 指数 
标准 化 .无量 纲 的 分 布 指数 "站 可 消除 不 同 地 形 
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图 2 不 同窗 口 面积 下 5-5 变化 范围 散 点 图 


Fig.2 Scatter plot of S—S; variation range under different 


window areas 


子 分 级 .山脉 冰川 面积 差异 的 影响 ,能 准确 客观 
揭示 青藏 高 原 不 同 山 脉 冰川 在 各 地 形 因 子 上 的 实 
际 分 布 情况 与 标准 分 布 的 偏差 ,其 计算 公式 如 下 : 
P=(S,/S) x(S/S,) (3) 

式 中 :P 代 表 分 布 指 数 ;e 代 表 地 形 因 子 ;5 代表 在 e 
地 形 因 子 特定 等 级 下 的 i 座 山 冰川 面积 ;5; 代 表 i 座 
山 冰川 总 面积 ;5 为 整个 研究 区 的 面积 ;5. 为 整个 区 
域 e 地 形 因子 特定 等 级 下 的 总 面积 。 分 布 指数 的 大 
小 对 于 冰川 在 地 形 上 的 空间 分 布 有 着 一 定 的 指示 
意义 :P>1 时 ,表示 冰川 的 分 布 呈现 优势 分 布 ,说 明 
该 地 形 带 的 条 件 适宜 冰川 的 生存 或 发 育 ;P<1 时 , 呈 
现 劣势 分 布 ,说 明 该 地 形 级 别 不 适宜 冰川 生存 或 发 
育 。 优 势 度 随 着 P 值 增 大 而 提高 ,反之 则 劣势 度 加 
剧 。 男 外 , 千 冰 川 分 布 曲线 比较 平缓 ,起 伏 不 大 , 表 
明 冰 川 的 分 布 与 标准 分 布 偏离 较 小 ,对 地 形 差异 的 
适宜 性 较 广 ;车 曲线 起 伏 较 大 , 则 表明 冰川 对 地 
具有 较 强 选择 性 ,在 其 优势 地 形 分 级 上 发 育 较 多 。 
2.3 地 理 探 测 器 

地 理 探 测 需 是 探测 空间 分 异性 ,以 及 揭示 其 背 
后 驱动 因子 的 一 种 新 的 统计 学 方法 ,其 因子 探测 gq 
统计 量 可 以 度量 空间 分 异性 \ 探 测 解 释 因 子 中 ,其 公 
式 如 下 : 


NS 


L 3 
2 No; 
h=1 


Ime (4) 
AXUPihcl, LAA Seek AE aNE ; NAN 2; 
别 为 层 h 和 全 区 的 单元 数 ; o, 和 为 自 变量 驱动 因 
子 和 因 变 量 冰川 分 布 的 方差 。g 的 值 域 为 [0,1], 值 越 
大 表示 驱动 因子 对 冰川 分 布 的 解释 力 越 强 ,反之 则 越 
弱 。 极 端 情况 下 ,g 值 为 1 表明 驱动 因子 完全 控制 了 
冰川 空间 分 布 ,g 值 为 0 表明 二 者 之 间 没 有 任何 关系 。 


3 结果 与 分 析 


3.1 基于 地 形 梯度 的 青藏 高 原 冰川 分 布 格局 

为 探究 冰川 在 各 地 形 因 子 上 的 分 布 特征 ,将 青 
藏 高 原 地 形 起 伏 度 (0.3989~11.5973) 海拔 (297-~ 
8781 m) 和 坡度 (0°~86.3963° ) 按 自然 断 点 法 分 为 9 
级 ,将 坡 向 按 无 坡 向 北 坡 、 东 北 坡 、 东 坡 、 东 南 坡 、 
南 坡 、 西 南 坡 、 西 坡 、 西 北 坡 分 为 9 级 ,然后 分 别 计算 
青藏 高 原 冰 川 在 地 形 起 伏 度 海拔 . 坡 向 .坡度 的 分 
布 指 数 (图 3)。 由 图 3 可 知 ,青藏 高 原 冰 川 在 地 形 起 
IRE 海拔 上 的 分 布 曲 线 起 伏 较 大 , 且 随 着 地 形 起 
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图 3 各 地 形 因子 冰川 的 分 布 特征 


Fig. 3 The distribution characteristics of glaciers with different topographic factors 


伏 度 ,海拔 的 升 高 ,青藏 高 原 冰 川 的 分 布 指数 呈现 
递增 的 趋势 ,并 在 地 形 起 伏 度 第 6-9 分 级 (5.0795~ 
11.5973) 海拔 第 8~9 分 级 (5082~8781 m) EE MÈ 
势 分 布 ,说 明 冰 川 对 地 形 起 伏 度 ,海拔 具有 较 强 选 
择 性 ; 青藏 高 原 冰川 在 坡 向 .坡度 上 的 分 布 曲 线 比 
较 平 缓 ,起 伏 不 大 ,表明 冰川 对 坡 向 .坡度 的 适宜 性 
较 广 。 冰 川 在 北 坡 .东北 坡 . 东 坡 和 西北 坡 上 呈现 
优势 分 布 ,其 中 , 北 坡 分 布 指数 最 高 ,其 次 是 东北 
JE ;另外 ,冰川 在 坡度 第 2~7 分 级 (5.0623°~ 
36.1110°) 上 的 曲线 特别 平缓 ,集中 在 标准 1 的 附近 、 
在 7~9 分 级 (36.1110°~86.3963°) 上 , 随 着 坡度 变 陡 ， 


分 布 指数 呈现 递增 的 趋势 。 
3.2 地 形 对 冰川 分 布 影响 的 定量 分 析 

地 理 探测 器 因子 探测 可 以 明确 地 形 因 子 对 冰 
川 分 布 空间 分 异 的 影响 程度 。 利 用 ArcGIS 软 件 , 生 
成 覆盖 青藏 高 原 全 区 的 2 kmx2 km 的 采样 网 格 ,并 
选取 距离 青藏 高 原 冰 川 5 km 范围 内 的 网 格 ,共计 
179736 个 。 然 后 ,计算 每 个 网 格 冰川 的 面积 分 布 百 
分 比 ,以 其 作为 因 变 量 Y, 并 将 其 赋 给 网 格 的 中 心 
点 。 将 海拔 .地 形 起 伏 度 .坡度 . 坡 向 4 个 地 形 因 子 
作为 自 变量 X, 借 助 地 理 探测 器 技术 计算 的 因子 探 
测 及 交互 探测 结果 见 表 1。 


表 1 地 形 因 子 对 冰川 分 布 的 影响 


Tab.1 Influence oftopographic factors on glacier distribution 


妹子 探测 g 值 交互 探测 
w OE guy apu mnene ange PANIE EE Waag 

阿尔 金山 0.5901 0.3296 0.0158 0.0172 0.6651 0.6258 0.1082 0.6221 0.4106 0.4834 
祁连山 0.5205 0.2006 0.0046 0.0123 0.5480 | 0.5359 0.0343 0.5596 0.2277 0.2792 
青 古 拉 山 0.4717 0.2658 0.0017 0.0028 0.4794 0.5236 0.0162 0.6128 0.2770 0.4488 
BELL 0.4346 — 0.2242 0.0128 0.0070 0.4498 — 0.4756 | 0.0252 0.4610 0.2397 0.4007 
兴 都 库 什 山 0.3885 0.2045 0.0454 0.0108 0.4077 . 04451 0.0657 0.4326 0.2216 0.2951 
自 米尔 0.3445 0.2319 0.0279 0.0105 0.3632 0.3920 0.0448 0.3815 0.2496 0.3109 
横断 山 0.2400 0.1328 0.0152 0.0056 0.2507 0.2631 0.0342 0.2829 0.1511 0.1865 
羌 塘 高 原 0.1851 0.2220 0.0003 0.0044 0.1920 0.2568 0.0237 0.2428 0.2325 0.3218 
喜马拉雅 山 0.1720 0.1256 0.0422 0.0039 0.1768 0.2115 0.0485 0.2150 0.1330 0.1979 
喀 喇 昆仑 山 0.1593 0.1508 0.0475 0.0053 0.1668 0.2475 0.0564 0.2493 0.1612 0.2520 
念 青 唐古拉 山 0.1201 0.0915 0.0202 0.0028 0.1258 0.1556 0.0290 0.1733 0.0989 0.1357 
冈 底 斯 山 0.0673 0.0992 0.0119 0.0019 0.0698 0.0811 0.0286 0.1170 0.1082 0.1213 
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由 表 1 可 知 ,各 山脉 冰川 的 海拔 、 地 形 起 伏 度 的 
4 值 远 高 于 坡度 和 坡 向 。 且 除尘 塘 高 原 .四 底 斯 山 
外 ,其 余 山 脉 冰川 的 g 值 均 呈 现 出 海拔 > 地 形 起 伏 度 
的 特征 。 男 外 ,坡度 和 坡 向 的 g 值 均 很 小 ,说 明 坡度 
和 坡 向 对 冰川 空间 分 布 的 影响 力 明 显 偏 弱 。 可 见 ， 
海拔 和 地 形 起 伏 度 是 控制 青藏 高 原 各 山脉 冰川 空 
间 分 布 的 主导 因素 。 

交互 探测 结果 表明 ,各 地 形 因子 对 冰川 分 布 格 
局 的 影响 并 不 独立 ,任意 两 个 地 形 因 子 交 互 作用 均 
高 于 单 因 子 作 用 , 均 呈 现 双 因 子 增强 和 非 线 性 增 
强 , 不 存在 独立 及 减弱 的 关系 。 这 也 说 明 本 文选 取 
的 4 个 影响 因子 中 任意 两 个 因子 的 交互 作用 对 青藏 
高 原 冰 川 分 布 格局 的 影响 更 为 显著 。 海 拔 和 地 形 
起 伏 度 的 共同 作用 对 冰川 分 布 的 影响 最 大 。 虽 然 ， 
坡度 和 坡 向 的 单 因子 q 值 都 较 低 ,但 与 其 他 因子 交 
互 后 均 超过 单 因子 单独 作用 对 冰川 分 布 格局 的 影 
响 , 因 此 ,这 些 因 子 的 影响 也 不 可 忽视 。 


另外 ,同一 地 形 因子 或 地 形 因子 组 合 对 各 山脉 
冰川 分 布 空间 分 异 的 解释 力 强 弱 存 在 明显 的 差 
异 。 其 中 ,海拔 对 阿尔 金山 、 祁 连 山 、 唐 古 拉 山 \ 昆 
仑 山 、 兴 都 库 什 山 等 冰川 空间 格局 分 异 的 解释 力 较 
强 ,g 值 均 超过 0.3; 而 对 横断 山 、 羌 塘 高 原 、 喜 马 拉 雅 
山 、 喀 喇 昆 仑 山 、 念 青 唐古拉 山 等 冰川 空间 格局 分 
异 的 解释 力 较 弱 ,对 冈 底 斯 山 冰 川 空间 格局 分 异 的 
解释 力 最 弱 ,9 值 仅 为 0.0673 ,上 且 冈 底 斯 山 的 冰川 分 
布 受 各 地 形 因 子 及 其 组 合 的 影响 均 较 小 ,数值 均 小 
3-0.12 
3.3 地 形 起 伏 度 海拔 梯度 上 各 山脉 冰川 分 布 特征 

为 了 进一步 探究 青藏 高 原 各 山脉 冰川 在 地 形 
起 伏 度 ,海拔 梯度 上 的 分 布 特征 ,将 地 形 起 伏 度 海 
拨 按 等 距 法 分 为 50 级 ,运用 分 布 指数 ,分别 计算 青 
藏 高 原 各 山脉 冰川 在 地 形 起 伏 度 海拔 梯度 上 的 分 
布 指数 (图 4, 图 5)。 

由 图 4 可 知 ,各 山脉 冰川 在 地 形 起 伏 度 梯度 上 
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图 4 各 山脉 冰川 在 地 形 起 伏 度 梯度 上 的 分 布 


Fig. 4 Distribution map of mountain glaciers based on RDLS 
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图 5 各 山脉 冰川 在 海拔 梯度 上 的 分 布 


Fig. 5 The distribution of glaciers in each mountain range on the elevation gradient 


BRUSRWIEAS. BSH ERE EEG LU .由 
米尔 等 冰川 优势 分 布 范 于 在 20 个 区 间 以 上 , 且 在 优 
势 分 布 区 间 分 布 指数 值 整体 较 高 。 其 中 , 喜 马 拉 和 雅 
山脉 冰川 优势 分 布 范 围 最 广 ,为 23~49 区 间 
(5.3262~11.3733) , 且 在 该 范围 分 布 指数 随 着 地 形 
起 伏 度 增加 总 体 呈 增长 趋势 。 喀 喇 昆仑 山脉 冰川 
优势 分 布 范围 次 之 ,为 23-48 IX [8] (5.3262~ 
11.1494) ,上 且 随 着 地 形 起 伏 度 的 增加 ,分 布 指数 呈现 
先 增长 后 又 降 的 趋势 。 兴 都 库 什 山 、 横 断 山 . 念 青 
唐古拉 出 的 冰川 优势 分 布 区 间 在 17~20 个 区 间 ,分 
布 指数 随地 形 梯度 的 形态 呈现 扁平 型 ;昆仑 山 .内 
底 斯 山 、 唐 古 拉 山 、 阿 尔 金 山 、 祁 连 山 、 羌 塘 高 原 的 
冰川 优势 分 布 区 间 在 10 个 区 间 以 内 , 均 旦 现 瘦 高 
型 , 且 多 为 双 峰 、 多 峰 的 形态 ,其 中 , 羌 塘 高 原 冰川 
优势 分 布 区 间 最 罕 。 


由 图 5 可 知 , 各 山脉 冰川 面积 在 海拔 梯度 上 的 
分 布 格局 差异 显著 ,上 且 形 态 各 异 。 喜 马 拉 雅 山脉 冰川 
优势 分 布 范围 最 广 , 为 29~50 区 间 (5048~8781m) , H. 
随 着 海拔 的 升 高 ,呈现 急剧 增加 ,后 趋 于 平稳 的 趋 
势 。 喀 喇 昆 仑 山脉 冰川 优势 分 布 区 间 较 广 ,为 29~ 
49 X [B] (5048~8611 m) ,在 30~40 区 间 呈 现 剧 增 的 趋 
势 ,随后 呈现 递减 趋势 。 其 余 10 座 山脉 冰川 的 优势 
分 布 均 呈 近似 正 态 分 布 。 兴 都 库 什 出 、 帕 米尔 高 原 
和 横断 山 ,优势 分 布 区 间 分 别 为 23~44 区 间 (4369~ 
7763 m) ,277-43 区 间 (4709~7593 m) Fil 29~43 [X [R] 
(5048~7593 m); 念 青 唐 古 拉 山 .昆仑 山 ERTL 
唐古拉 山 ` 阿 尔 金山 . 羌 塘 高 原 .祁连山 的 冰川 优势 
分 布 区 间 数 目 在 15 个 区 间 以 内 ,分 布 指数 曲线 呈现 
瘦 高 型 正 态 分 布 。 其 中 ,祁连山 冰川 优势 分 布 区 间 
最 罕 ,优势 分 布 区 间 数 目 仅 为 5 个 区 间 。 
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-— 局 存在 明显 差异 ,峰值 的 高 度 不 同 , 且 优势 分 布 的 
4 讨论 范围 也 不 同 ,进而 呈现 出 各 异 的 形态 特征 。 较 低 纬 


青藏 高 原 冰川 的 分 布 格局 与 地 形 因 子 之 间 的 
关系 一 直 是 地 理学 和 冰川 学 研究 的 重要 议题 。 本 
研究 发 现 ,青藏 高 原 冰川 对 坡 向 坡度 的 适宜 性 较 
广 ,对 地 形 起 伏 度 海拔 的 选择 性 较 强 , 且 在 高 地 形 
起 伏 度 、 高 海拔 上 呈现 优势 分 布 。 地 理 探测 器 因子 
探测 结果 也 表明 ,在 所 选 的 4 个 地 形 因子 中 ,海拔 和 
地 形 起 伏 度 是 控制 冰川 分 布 的 主要 因素 , 且 交 互 探 
测 结果 表明 二 者 共同 作用 对 冰川 空间 分 异性 影响 
最 大 。 这 与 前 人 的 研究 结论 ““ 相 一 致 。 但 同一 地 
形 因 子 或 地 形 因子 组 合 对 各 山脉 冰川 分 布 空间 分 
异 的 解释 力 强 弱 存 在 明显 的 差异 ,这 可 能 与 各 山脉 
独特 的 区 域 位 置 与 气候 条 件 有 关 。 青 藏 高 原 冰 川 
分 布 受到 多 种 气象 因素 的 影响 ,其 中 温度 和 降水 是 
决定 冰川 发 育 的 物质 基础 ”。 这 些 气象 因素 与 地 
形 因子 共同 作用 ,决定 了 冰川 的 形成 .发 育 和 退缩 。 
青藏 高 原 冰 川 主要 受 西南 季风 、 东 南 季风 南亚 季风 等 
季风 环流 和 西风 环流 及 其 相互 作用 的 影响 "。 各 山脉 
外 侧 迎 风 坡 降水 丰富 ,有 利于 冰川 发 育 ,又 因 高 大 山体 
遮挡 作用 ,背风 坡 降水 较 少 。 青 藏 高 原 亚洲 夏季 风气 
候 态 北 界 沿 着 祁连山 东 段 -唐古拉 山 - 羌 塘 高 原 - 内 底 
斯 山 - 西 喜马拉雅 山 分 布 ”。 该 界线 以 北 为 西风 影响 
区 ,降水 稀少 ,气候 干 冷 , 地 形 因 子 对 该 区 较 多 山脉 
冰川 分 布 空间 分 异 的 解释 力 较 强 ;该 界线 以 南 为 季 
影响 区 ,降水 较 多 ,气候 湿润 ,地 形 因子 对 该 区 较 
多 山脉 冰川 分 布 空 间 分 异 的 解释 力 较 弱 ; 而 受 两 个 
环流 系统 共同 影响 的 西风 -季风 过 渡 区 则 更 为 复 
杂 。 其 中 ,由 底 斯 山 - 念 青 唐古拉 山 又 为 青藏 高 原 
年 均 温 0% 和 最 冷 月 均 温 -10 的 分 界线 ,其 南北 两 
侧 的 水 热 状况 差异 显著 , 南 侧 气候 较 北 部 温暖 湿 
润 , 且 温 度 差 异 呈 现 * 跃 变 " 特 征 ” ,从 而 决定 了 该 
界线 两 侧 冰川 分 布 格局 的 差异 。 其 独特 的 气候 特 
征 可 能 是 导致 地 形 因 子 对 内 底 斯 山 - 念 至 唐古拉 山 
冰川 分 布 空间 分 异 的 解释 力 最 弱 的 原因 。 

地 形 和 气候 共同 作用 是 导致 不 同 山脉 冰川 在 
地 形 梯 度 上 分 布 格局 差异 的 原因 。 此 外 ,青藏 高 原 
各 山脉 冰川 在 地 形 梯 度 上 的 分 布 多 呈现 近似 正 态 
分 布 特征 , 即 在 某 一 个 地 形 带 上 分 布 指数 出 现 峰 
值 ,表征 该 地 形 带 的 气候 或 地 形 最 适宜 冰川 的 生存 
和 发 展 。 而 不 同 山脉 冰川 在 地 形 梯度 上 的 分 布 格 


度 地 区 的 山脉 冰川 ,如 喜马拉雅 山 ` 念 青 唐古拉 山 、 
喀 喇 昆仑 山 , 因 所 处 区 域 热量 充足 ,温度 高 的 原因 , 
冰川 只 能 在 海拔 较 高 的 范围 存在 ,另外 ,降水 相对 
BEE ,保障 了 冰川 向 较 低 海拔 延伸 , 故 冰 川 优势 分 
布 区 间 较 广 ;内 底 斯 山 ORE RR .昆仑 山 的 冰川 类 
型 多 为 极 大 陆 型 冰川 ,所 处 区 域 降水 稀少 ,冰川 赖 
以 生存 的 物质 条 件 相 对 缺乏 ,冰川 延展 范围 较 小 ， 
故 分 布 海拔 相对 集中 ,分布 特 征 呈 现 瘦 高 型 ”…。 另 
外 , 因 其 冰川 厚度 大 , 受 地 形 起 伏 的 影响 相对 较 弱 ， 
导致 呈现 多 峰 的 状态 ; 念 青 唐古拉 山东 段 横断 山 
脉 冰 川 属于 海洋 型 冰川 ,地 处 中 国 第 一 、 二 阶梯 的 
过 渡 带 ,位 于 东南 和 西南 两 大 季风 交汇 区 , 受 区 域 
大 气 环流 及 地 形 因素 综合 影响 ,降水 充沛 ,冰川 得 
以 向 较 低 海拔 延伸 , 故 冰川 优势 分 布 区 间 较 广 。 又 
因 地 形 复杂 地域 差异 性 强 ,导致 呈现 多 峰 的 状态 。 

综 上 所 述 ,地 形 因 子 对 冰川 分 布 的 指示 能 力 是 
显而易见 的 ,但 在 某 些 特 定 的 区 域 ,这 种 指示 能 力 
可 能 会 受到 气温 、 降 水 、 风 向 、 区 位 等 因素 的 影响 。 
本 文 仅 探 究 了 青藏 高 原 冰川 面积 在 各 地 形 梯度 上 
的 分 布 格局 ,为 了 深入 人 研究 基于 地 形 梯 度 的 冰川 分 
布 格局 ,还 应 开展 对 冰川 体积 .厚度 .长 度 .平衡 线 
等 方面 的 研究 。 又 因 人 缺乏 青藏 高 原 地 区 气象 观测 
数据 ,对 青藏 高 原 冰 川 影响 成 因 的 分 析 沿 不 充分 。 
在 未 来 的 研究 中 , 需 而 合 不 同 尺 度 冰川 的 分 布 格局 
和 时 空 变化 规律 ,综合 考虑 气象 地形、 人文 等 多 种 
影响 因素 ,建立 冰川 系统 模型 ,以 模拟 预测 冰川 的 
变化 趋势 ,并 定量 评估 冰川 变化 对 全 球 气候 与 环境 
变化 的 影响 。 


5 结论 


本 文 基于 第 6 版 Randolph 冰川 编目 和 NASA- 
DEM 数据 ,采用 均值 变 点 法 确定 了 青藏 高 原 地 区 地 
形 起 伏 度 的 最 佳 计算 窗口 ,并 利用 分 布 指数 和 地 理 
探测 器 技术 ,分 析 了 青藏 高 原 各 山脉 冰川 在 地 形 起 
伏 度 海拔 , 坡 向 .坡度 上 的 分 布 格局 及 其 成 因 。 主 
要 结论 如 下 : 

CL) 青藏 高 原 冰川 对 坡 向 .坡度 的 适宜 性 较 广 ， 
对 地 形 起 伏 度 ,海拔 的 选择 性 较 强 ,上 且 在 高 地 形 起 
伏 度 .高 海拔 上 呈现 优势 分 布 。 随 着 地 形 起 伏 度 、 


海拔 的 升 高 ,青藏 高 原 冰 川 的 分 布 指数 呈现 递增 的 
趋势 , 且 在 地 形 起 伏 度 第 6-9 海拔 第 8-9 分 级 上 呈 
现 优 势 分 布 ;青藏 高 原 冰川 在 北 坡 、 东 北 坡 、 东 坡 和 
西北 坡 上 呈现 优势 分 布 , 其 中 , 北 坡 分 布 指 数 最 高 ， 
其 次 是 东北 坡 ; 青 藏 高 原 冰川 对 坡度 的 适宜 性 最 
广 , 在 坡度 第 2~7 分 级 分 布 指数 集中 在 标准 1 的 附 
近 , 在 7~9 分 级 上 , 随 着 坡度 变 陡 , 分 布 指数 呈现 递 
增 的 趋势 。 

(2) 地 形 对 冰川 的 发 育 影响 较 大 ,各 地 形 因 子 
对 冰川 空间 分 布 的 影响 显著 不 同 ,海拔 和 地 形 起 伏 
度 是 控制 冰川 分 布 的 主要 因素 ,其 次 是 坡度 和 坡 
向 。 交 互 探测 结果 表明 ,海拔 坡度 . 坡 向 和 地 形 起 
伏 度 对 冰川 分 布 格局 的 影响 并 不 独立 ,任意 两 个 地 
形 因子 对 冰川 分 布 的 交互 作用 高 于 单 因子 作用 , 均 
呈现 双 因 子 增 强 和 非 线 性 增强 ,不 存在 独立 及 减弱 
的 关系 。 海 拔 和 地 形 起 伏 度 的 共同 作用 对 冰川 空 
间 分 异性 影响 最 大 。 

(3) 在 地 形 起 伏 度 分 级 上 ,喜马拉雅 山 . 喀 喇 昆 
仑 山 .由 米尔 高 原 的 冰川 优势 分 布 范 围 在 20 个 区 间 
以 上 , 且 在 优势 分 布 区 间 分 布 指数 值 整体 较 高 。 其 
中 ,喜马拉雅 山脉 冰川 优势 分 布 范 围 最 广 。 在 海拔 
分 级 上 ,喜马拉雅 山脉 冰川 优势 分 布 范围 最 广 , 且 
随 着 海拔 升 高 呈现 急剧 增加 ,后 趋 于 平稳 的 趋势 。 
喀 喇 昆仑 山脉 冰川 优势 分 布 区 间 较 广 ,呈现 先 剧 增 
后 递减 的 趋势 。 其 余 10 座 山脉 冰川 的 优势 分 布 均 
呈 近 似 正 态 分 布 , 但 不 同 山脉 的 形态 差异 显著 。 
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Distribution pattern and causes of glaciers in the Tibetan Plateau based 


on terrain gradient 
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475004, Henan, China) 


Abstract: The formation and development of glaciers depend on a combination of topographic conditions and 
climatic factors, such as precipitation and temperature. However, the mechanisms underlying glacier distribution 
and variation remain uncertain. Our study aims to explore the relationship between glacier distribution patterns 
and topographical factors and uncover the driving forces behind them. To accomplish this, we utilized the 
Randolph Glacier Inventory version 6 and NASADEM for our analysis. Initially, we derived four key 
topographical factors based on NASADEM: relief degree of land surface (RDLS), altitude, slope, and aspect. 
Because the calculation of RDLS has scale dependence, the key is to determine the optimal analysis window 
using the average change- point analysis approach. The identified optimal analysis window size for RDLS 
evaluation in the study area was a rectangular neighborhood of 36 x 36 pixels, corresponding to an area of 
approximately 1.17 km’. Subsequently, we examined glacier distribution patterns across various terrains in the 
Tibetan Plateau using the distribution index. This index elucidates the disparities between actual and standard 
glacier distributions across various mountain chains, excluding area disturbances. Finally, we employed the 
geodetector method to quantitatively assess the spatial interplay between glacier distribution patterns and 
topographical elements. This innovative statistical approach identifies spatially stratified heterogeneity, pinpoints 
explanatory factors, and evaluates interactive relationships between variables. The results show that: (1) Glacier 
distribution in the Tibetan Plateau shows strong selectivity for RDLS and altitude, with wide suitability for aspect 
and slope. The predominant distribution involves large RDLS and high altitudes. (2) Topography significantly 
influences glacier development, with varied effects of topographic factors on the spatial distribution of glaciers. 
Altitude and RDLS emerged as dominant factors controlling the distribution of glaciers, followed by slope and 
aspect. The interactive detection revealed that the combined effect of altitude and RDLS had a dominant impact 
on spatially stratified glacier heterogeneity. (3) Among altitude and RDLS gradients, the Himalayas Mountains 
had the most extensive glacier distribution, followed by the Karakoram Mountains. The other ten mountain 
ranges demonstrated diverse glacier distributions, though they mostly adhered to a normal pattern. This study will 
serve as a theoretical reference and provide data support for regional hydrology research. It offers vital guidance 
for appropriate water resource management in arid lands. 


Keywords: glacier distribution; terrain gradient; distribution index; geodetector; Tibetan Plateau 


